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Abstract. A major problem in Cosmic Microwave Back- 
ground (CMB) anisotropy mapping, especially in a total- 
power mode, is the presence of low-frequency noise in the 
data streams. If unproperly processed, such low-frequency 
noise leads to striping in the maps. To deal with this prob- 
lem, solutions have already been found for mapping the 
CMB temperature fluctuations but no solution has yet 
been proposed for the measurement of CMB polarization. 
Complications arise due to the scan-dependent orientation 
of the measured polarization. In this paper, we investigate 
a method for building temperature and polarization maps 
with minimal striping effects in the case of a circular scan- 
ning strategy mission such as the Planck mission. 

Key words: methods: data analysis; cosmology: cosmic 
microwave background; polarization 
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K*" ■ 1. Introduction 



Theoretical studies of the CMB have shown that the ac- 
curate measurement of the CMB anisotropy spectrum Cf 
with future space missions such as Planck will allow 
for tests of cosmological scenarios and the determination 
of cosmological parameters with unprecedented accuracy. 
Nevertheless, some near degeneracies between sets of cos- 
mological parameters yield very similar CMB temperature 
anisotropy spectra. The measurement of the CMB polar- 



ization and the computation of its power spectrum (Seljak 



1996; Zaldarriaga, 1998) may lift to some extent some of 



these degeneracies. It will also provide additional infor- 
mation on the reionization epoch and on the presence of 
tensor perturbations, and may also help in the identifi- 
cation and removal of polarized astrophysical foregrounds 



19991 ) 



( Kinney, 1998 ; Kamionkowski, 1998| ; Prunet and Lazarian, 



A successful measurement of the CMB polarization 
stands as an observational challenge; the expected polar- 
ization level is of the order of 10% of the level of tem- 
perature fluctuations (AT/T ~ 10~^). Efforts have thus 



gone into developing techniques to reduce or eliminate 
spurious non-astronomical signals and instrumental noise 
which could otherwise easily wipe out real polarization sig- 



nals. In a previous paper (Couchot et al., 1999), we have 



shown how to configure the polarimeters in the focal plane 
in order to minimize the errors on the measurement of the 
Stokes parameters. In this paper, we address the problem 
of low frequency noise. 

Low frequency noise in the data streams can arise due 
to a wide range of physical processes connected to the 
detection of radiation, l/f noise in the electronics, gain 
instabilities, and temperature fluctuations of instrument 
parts radiatively coupled to the detectors, all produce low 
frequency drifts of the detector outputs. The spectrum of 
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Fig. 1. The power spectrum of the K34 bolometer from 
Caltech, the same type of bolometer planned to be used 
on the Planck mission. The measurement was performed 
at the SYMBOL test bench at I.A.S., Orsay (supplied by 
Michel Piat). The knee frequency of this spectrum is ~ 
0.014 Hz and the planned spin frequency for Planck is 
0.016 Hz. We can model the spectrum as the function 
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the total noise can be modeled as a superposition of white 
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noise and components behaving like 1//" where a> 1, as 
shown in Fig. |^. 

This noise generates stripes after reprojection on maps, 
whose exact form depends on the scanning strategy. If not 
properly subtracted, the effect of such stripes is to degrade 
considerably the sensitivity of an experiment. The elimina- 
tion of this "striping" may be achieved using redundancies 
in the measurement, which are essentially of two types for 
the case of Planck: 

— each individual detector's field of view scans the sky 
on large circles, each of which is covered consecutively 
many times 60) at a rate of about /spin 1 rpm. 
This permits a filtering out of non scan-synchronous 
fluctuations in the circle constructed from averaging 
the consecutive scans. 

— a survey of the whole sky (or a part of it) involves 
many such circles that intersect each other (see Fig. 
^); the exact number of intersections depends on the 
scanning strategy but is of the order of 10^ for the 
Planck mission: this will allow to constrain the noise 
at the intersection points. 




Fig. 2. The MoUweide projection of 3 intersecting circles. 
For clarity, the scan angle between the spin axis and the 
main beam axis is set to 60° for this figure. 



One of us ( Delabrouille, 199^ ) has proposed to remove 
low frequency drifts for unpolarized data in the frame- 
work of the Planck mission by requiring that all mea- 
surements of a single point, from all the circles intersect- 
ing that point, share a common sky temperature signal. 
The problem is more complicated in the case of polarized 
measurements since the orientation of a polarimeter with 
respect to the sky depends on the scanning circle. Thus, a 
given polarimeter crossing a given point in the sky along 
two different circles will not measure the same signal, as 
illustrated in Fig. ||. 

The rest of the paper is organized as follows: in Sect. ||, 
we explain how we model the noise and how low frequency 
drifts transform into offsets when considering the circles 
instead of individual scans. In Sect. we explain how 
polarization is measured. The details of the algorithm for 




Fig. 3. The orientation of polarimeters at an intersection 
point. This point is seen by two different circles corre- 
sponding to two different orientations of the polarimeters 
in the focal plane. For clarity, we have just represented 
one polarimeter. 

removing low- frequency drifts are given in Sect. ^. We 
present the results of our simulations in Sect. ^ and give 
our conclusions in Sect. ^. 

2. Averaging noise to offsets on circles 

As shown in Fig. |l|, the typical noise spectrum expected 
for the Planck High Frequency Instrument (HFI) fea- 
tures a drastic increase of noise power at low frequencies 
/ < 0.01 Hz. We model this noise spectrum as: 
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(1) 



The knee frequency /knee is defined as the frequency 
at which the power spectrum due to low frequency contri- 
butions equals that of the white noise. The noise behaves 
as pure white noise with variance at high frequencies. 
The spectral index of each component of the low-frequency 
noise, a^, is typically between 1 and 2, depending on the 
physical process generating the noise. 

The Fourier spectrum of the noise on the circle ob- 
tained by combining N consecutive scans depends on the 
exact method used. The simplest method, setting the cir- 
cle equal to the average of all its scans, efficiently filters 
out all frequencies save the harmonics of the spinning fre- 



quency (Delabrouille et al., 1998b). Since the noise power 
mainly resides at low frequencies (see Fig. |l|), the averag- 
ing transforms - to first order - low frequency drifts into 
constant offsets different for each circle and for each po- 
larimeter. This is illustrated in the comparison between 
Figs. |4| and ||. More sophisticated methods for recombin- 
ing the data streams into circles can be used, as min- 
imization, Wiener filtering, or any map-making method 
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projecting about 6 x 10^ samples onto a circle of about 
5 X 10'^ points. For simplicity, we will work in the follow- 
ing with the circles obtained by simple averaging of all its 
consecutive scans. 
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ing white noise of the h polarimeters is described by one 
constant h x h matrix. 

3. The measurement of sky polarization 

3.1. Observational method 

The measurement with one polarimeter of the linear po- 
larization of a wave coming from a direction n on the 
sky, requires at least three measurements with different 
polarimeter orientations. Since the Stokes parameters Q 
and U are not invariant under rotations, we define them 
at each point n with respect to a reference frame of tan- 
gential vectors {e\, ep). The output signal given by a po- 
larimeter looking at point h is: 

Mpoiar = i (/ -f Q COS 2* + C/ sin 2^-) (2) 

where VE* is the angle between the polarimeter and caQ. In 
the following, we choose the longitude-latitude reference 
frame as the fixed reference frame on the sky (see Fig. 0). 



Fig. 4. Typical 1//^ low frequency noise stream. Here, 
/knee = /spin = 0.016 Hz, a = 2 and a = 21 /xK (see Eq. 
1). This noise stream corresponds to 180 scans or 3 circles 
(60 scans per circle) or a duration of 3 hours. 
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Fig. 5. The residual noise on the 3 circles after averag- 
ing. To first approximation, low frequency drifts are trans- 
formed into offsets, different for each circle and each po- 
larimeter. Note the expanded scale on the y-axis as com- 
pared to that of Fig. ^. 

We thus model the effect of low frequency drifts as a 
constant offset for each polarimeter and each circle. This 
approximation is excellent for /knee < /spin- The remain- 




Fig. 6. The reference frame used to define the Stokes pa- 
rameters and angular position ^ of a polarimeter. "if lies 
in the plane {e\, e/3). 



3.2. Destriping method 

The destriping method consists in using redundancies at 
the intersections between circle pairs to estimate, for each 
circle i and each polarimeter p, the offsets Of on polarime- 
ter measurements. For each circle intersection, we require 
that all three Stokes parameters in a fixed reference frame 
in that direction of the sky, as measured on each of the 
intersecting circles, be the same. A minimization leads 

^ We do not consider the V Stokes parameter since no net 
circular polarization is expected through Thomson scattering. 
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to a linear system whose solution gives the offsets. By 
subtracting these offsets, we can recover the Stokes pa- 
rameters corrected for low-frequency noise. 

3.3. Formalism 

We consider a mission involving n circles. The set of all cir- 
cles that intercept circle i is denoted by T{i) and contains 
Nx{i) circles. For any pair of circles i and j, we denote the 
two points where these two circles intersect (if any) by 
{«, J, (5}. In this notation i is the circle currently scanned, 
j the intersecting circle in set and 5 indexes the two 
intersections {5 = 1(— 1) indexes the first (second) point 
encountered from the northernmost point on the circle) 
so that the points {j,i,—S} and {i,j,S} on the sky are 
identical. 

The Stokes parameters at point {i,j,S}, with respect 
to a fixed global reference system, are denoted by a 
3— vector Sij^s, with 



S i.i. — S 




(3) 



At intersection {i,j,6}, the set of measurements by 
h polarimeters travelling along the scanning circle i is a 
/i— vector denoted by Mij^s, and is related to the Stokes 
parameters at this point by (see Eq. 0): 



where is the h x 3 matrix: 



(4) 



1 cos2^p{i,j,S) sin2*p(i,j,(5) 



yi cos2*,,(i,j,(S) sm2-^h{i,j,5) J 



^p{i,j,5) £ [0, tt] is the angle between the orientation 
of polarimeter p and the reference axis in the fixed global 
reference frame (see Fig. The matrix A^j.s can be fac- 
torised as 



The rotation matrix Rij.s brings the focal plane to its po- 
sition at intersection (5} when scanning along circle 
i: 

1 

cos 2*1 (i, J, (5) sin24'i(i,j,,5) ). (7) 
-sin2*i(i,j, (5) cos2*i(i, j, (5) 

4. The algorithm 

J^.l. The general case 

To extract the offsets from the measurements, we use a 
minimization. This relates the measurements Mij^s to 
the offsets Oi and the Stokes parameters S'ij.a, using the 
redundancy condition (^). In order to take into account 
the two contributions of the noise (see Sect. ^) and of the 
Stokes parameters (see Eq. ^) , we model the measurement 
as: 

Mijj = A.Ri.j^s Sij^s + Oi + white noise. (8) 
so that we write 

iVr^ {M.^j^s -0,~ AR^,J,s S,,jj) . (9) 

where Ni is the h x h matrix of noise correlation between 
the h polarimeters. 

Minimization with respect to Oi and Sij,s yields the 
following equations: 

N,-' ^^'^ '^R^,3.S 5»j,5) = 0, (10) 

jei{i),s=±i 

and 

R;Is A^N,-' {M.^j^s O^- AR^,J^s 5»j,5) + 
Rjl_sA^Nj-^ {Mj,,^^s - Oj ~ ARJ.^,^s S^.,,s) = 0. 

(11) 

We can work with a reduced set of transformed mea- 
surements and offsets which can be viewed as the Stokes 
parameters in the focal reference frame and the associated 
offsets which are the 3 dimensional vectors: 



•■Ai^j^s — '^^■R'i,j,s • 



(5) 



The constant /i x 3 matrix A characterizes the geometrical 
setup of the h polarimeters in the focal reference frame: 



/I 1 



A 



\ 



1 cos2Ap sin2Ap 



\ 1 cos 2 Aft sin2A,i / 



(6) 



where Ap is the angle between the orientations of po- 
larimeters p and 1, so we have = ^'i -f Ap and Ai = 0. 



§,.j- 5 = Xr A' N,- M.^j^s and 

where Xi^ Ni^'^ A. 

Eqs. (|l^ and ^ then simplify to: 

(Sij,<5 ~ ~ Ri,j,S Sij^s) — 0, 

J6l(i),5=±l 

and 



(12) 



(13) 



R 



(14) 
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Ri.j^s Sij^s in Eq. (|lj) can be solved for and the result 
inserted in Eq. (p^). After a few algebraic manipulations, 
one gets the following linear system for the offsets as 
functions of the data Sij^s'. 

iei(i),<5=±i 

[a, -H(z,j,(5) A, -S^,,,5 + ^(i,i,(5)§,,,,-5 =0, (15) 
where the rotation 



R{i,3,5) = R 



r: 



.-5 



(16) 



brings the focal reference frame from its position along 
scan j at intersection {i,j,S} to its position along scan 
i at the same intersection (remember that {i,j,S} — 
{j,i,-S}). Note that R{i,j,S) = Rij,i,-S)-^. In this 
linear system, we need to know the measurements of the 
polarimeters at the points {i,j,S} and {j,i,—6}. These 
two points on circles i and j respectively will unlikely cor- 
respond to a sample along these circles. So we have linearly 
interpolated the value of the intersection points from the 
values measured at sampled points. For a fixed circle i, 
this is a 3 X Nx(i) linear system. In Eq. ( |l5| ) , i runs from 1 
to n, therefore the total matrix to be inverted has dimen- 
sion 3n X 3n. However, because the rotation matrices are 
in fact two dimensional (see Eq. |^), the intensity compo- 
nents Af of the offsets only enter Eq. (jlj) through their 
differences Aj — so that the linear system is not invert- 
ible: the rank of the system is 3n — 1. In order to compute 
the offsets, we can fix the intensity offset on one particular 
scanning circle or add the additionnal constraint that the 
length of the solution vector is minimized. 

Once the offsets A^ are known, the Stokes parameters 
in the global reference frame [ex, e^j) at a generic sampling 
k of the circle i, labeled by {i, k} are estimated as 



Si^k — Rik i^iM ^ . 



(17) 



where Ri^k is the rotation matrix which transforms the 
focal frame Stokes parameters into those of the global ref- 
erence frame. 

The quantities 2>i^k are the Stokes parameters mea- 
sured in the focal frame of reference at this point and 
are simply given in terms of the measurements Mi^k (see 
Eq. H) by: 



X-^A^N-^M,^k- 



The matrix 



AT: 



Stokes 



(18) 



(19) 



is the variance matrix of the Stokes parameters on circle 
i. Note that this algorithm is totally independent of the 
pixelization chosen which only enters when reprojecting 
the Stokes parameters on the sphere. 



4-.2. Uncorrelated polarimeters, with identical noise 

When the polarimeters are uncorrelated with identical 
noise, the variance matrix reduces to Ni — \la^ and 
the matrices Xi can all be written as 



Xi = —X with X = A' A 



Case of "Optimized Configurations" 

We have shown ( Couchot et al., 1999| ) that the polarime- 
ters can be arranged in "Optimized Configurations" , 
where the h polarimeters are separated by angles of 7r//i. 
If the noise level of each of the h polarimeters is the same 
and if there are no correlation between detector noise, then 
the errors of the Stokes parameters are also decorrelated 
and the matrix X has the simple form: 



X = 




(20) 



Because this matrix commutes with all rotation matrices 
R{i,j,S), Eq. (O) simplifies further to 



N- 



1) A, 



E 



jGl(i),<5=±l ' 1 



§'i,i,<5 - -R(i, j, i5) §j 



(21) 



where the sum over b is explicit on the left side of the 
equation and we have defined an average error crx(i) along 
circle i by 

^l(i) _ \ - 2 

and where the rotation matrix R{i,j, S) is defined by Eq. 
(|l6|). Rotations R{i,j,S) and R{i,j,~6) — R{j,i,S)~^ 
correspond to the two intersections between circles i and 
j. Eq. ( |2l] ) can be simplified further. We can separate the 
Ai and the §ij,5 into scalar components related to the 
intensity: Af, j g and 2-vectors components related to 
the polarization: Af , §f^- g. We obtain then two separate 
equations, one for the intensity offsets Af, which is exactly 
the same as in the unpolarized case (see ^ : 



E 



A - A 



E ■ 



5), (22) 
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and one for the polarization offsets Aj 



N- 



E 



cos(2^',,) Af 



E ■ 

i6X(i),5=±l 



•j,i-S) ' 



(23) 



where = (5) — 5) is the angle of the 

rotation that brings the focal reference frame along scan 
j on the focal reference frame along scan i at intersection 
{i, j, S}. The two dimensional matrix TZ{i, j, S) is the rota- 
tion sub-matrix contained in R{i,j,6) (see Eq. ^6|). Note 
that all mixing between polarization components have dis- 
appeared from the left side of Eq. (p^)). Therefore we are 
left with two different n x n matrices to invert in order to 
solve for the offsets A, , instead of one Sn x 3n matrix. 



As in Eq. (|15| ), the linear system in Eq. ( |22D involves 
differences between the offsets Aj, the matrix is not in- 
vertible, and we find the solution in the same way as in 
the general case. On the other hand, for the polarized off- 
sets Af , the underlying matrix of Eq. (^3|) is regular as 
expected. 

5. Simulations and test of the algorithm 

5.1. Methods 

We now discuss how we test the method using numerical 
simulations. For each simulated mission, we produce sev- 
eral maps. The first, which we use as the standard "ref- 
erence" of comparison, is a projected map of a mission 
with only white noise, /knee = 0. The remaining maps in- 
clude 1// plus white noise streams with /knco — ?7/spin, 
T] G {1,2,5,10}. The first of these is an "untreated" im- 
age which is projected with no attempt to remove striping 
affects. In the second "zero-averaged" map, we attempt a 
crude destriping by subtracting its average to each circle. 
The final "destriped" map is constructed using the algo- 
rithm in this paper. We subtract the input maps (/, Q 
and U) from the final maps in order to get maps of noise 
residuals. Note that in case of a zero signal sky, setting 
the average of each circle to zero is better than destriping 
by nature because the offsets are only due to the noise. 
With a real sky, both signal and noise contribute to the 
average so that zeroing circle n ot onl y removes the noise 
but also the signal. Giard et al ( 1999 ) have attempted to 
refine their method by fitting templates of the dipole and 
of the galaxy before subtracting a baseline from each cir- 
cle. They concluded that an addit ional destriping (they 
used the algorithm of Delabrouille ( 1998| )) is needed. 



5.2. Simulated missions 

In order to test its efficiency, the destriping algorithm has 
been applied to raw data streams generated from simu- 
lated observations using various circular scanning strate- 



gies representative of a satellite mission as Planck, differ- 
ent "Optimized Configurations" , and various noise param- 
eters. The resulting maps were then compared with input 
and untreated maps to test the quality of the destriping. 

The input temperature (/) maps are the sum of 
galaxy, dipole, and a randomly generated standard CDM 
anisotropy map (we used HEALPIXQ and CMBfastQ). 
Similarly, the polarization maps Q and U are the sum 
of the galaxy and CMB polarizations. The CMB polariza- 
tion maps are randomly generated assuming a standard 
CDM scenario. For the galactic polarization maps, we con- 
structed a random, continuous and correlated vector field 
defined on the 2-sphere with a correlation length of 5° 
and a maximum polarization rate of 20% (100% gave sim- 
ilar results). Given the temperature map of the sky (not 
including CMB contribution), we can thus construct two 
polarization maps for Q and U . 

5.3. Results 

We first consider the case of destriping pure white noise 
and check that the destriping algorithm does not intro- 
duce spurious structure. Once this is verified, we apply 
the destriping algorithm to low frequency noise. We find 
that the quality of the destriping is significantly depen- 
dent on rj only. To demonstrate visually the quality of the 
destriping, we produce projected sky maps with the input 
galaxy, dipole and CMB signal subtracted. 

For temperature maps, we can compare Figs. 7,8, and 
9. The eye is not able to see any differences between the 
white noise map and the noise residual on the destriped 
map. We will see in the following how to quantify the pres- 
ence of structures. For the "zero-averaged circles" map, 
the level of the structure is very high and make it impos- 
sible to compute the power spectrum of the CMB (see Fig. 
13). 

For the Q Stokes parameter. Fig. 11 shows the de- 
striped map for 77 = 1. Fig. 12 shows a map where the off- 
sets are calculated as the average of each circle. The maps 
for U are very similar. As for the / maps, the destriped 
map is very similar to the white noise map. There ex- 
ist some residual structure on the "zero-averaged circles" 
map. To assess quantitatively the efficiency of the destrip- 
ing algorithm, we have first studied the power spectra Cj , 
Cf , Cf, Cj^ and Cj^ calculated from the /, Q and U 



1997) 



maps (Zaldarriaga and Seljak, 1997; Kamionkowski et al 



The reference sensitivity of our simulated mission is 
evaluated by computing the average spectra of 1000 maps 
of reprojected mission white noise. This reference sensitiv- 
ity falls, within sample variance, between the two dotted 
lines represented in Figs. 13, 14, 15 and 16 and below the 



littp : //www. tac .dk/'healpix/ 
tittp : //www. sns . ias . edu/~matiasz 
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Fig. 7. The MoUweide projection of the residuals of the 
/ — Stokes parameter for a white noise mission. The scale 
is in Kelvins. The parameters of the simulation leading to 
this map are described in 

Fig. 8. The MoUweide projection of the residuals of the 
/— Stokes parameter after destriping of 1// noise plus 
white noise with ?/ = 1 

Fig. 9. The MoUweide projection of the residuals of the 
/ — Stokes parameter after zeroing the average of the cir- 
cles, for 1// noise plus white noise with rj = 1. 

Fig. 10. The MoUweide projection of the residuals of the 
Q — Stokes parameter for a white noise mission. 

dotted line in Figs. 17 and 18. Figs. 13 and 14 show the 
spectra Cf corresponding to /knee/ /spin — 1 and 5 respec- 
tively, for the T field. Similarly, Figs. 15 and 16 are the 
spectra Cf. The B field is not represented because it is 
very similar to E. Figs. 17 and 18 represent the correlation 
between E and T: Cf^. For /knee//spin = 1, we see that 
we are able to remove very efficiently low frequency drifts 
in the noise stream: the destriped spectra obtained are 
compatible with the white spectrum (within sample vari- 
ance) . Similar quality destriping is achieved for any super- 
position of 1//, 1/ P and white noise, provided that the 
knee frequency is lower than or equal to the spin frequency. 
In the case of /kncc//spin = 5, the method as implemented 
here leaves some striping noise on the maps at low values 
of i. Modeling the noise as an offset is no longer ade- 
quate and a better model of the averaged low-frequency 
noise is required (superposition of sine and cosine func- 
tions for instance), or a more sophisticated method for 
constructing one circle from 60 scans. We again note that 
the value of the ratio /kncc//spin for both Planck HFI 
and LFI is likely to be very close to unity (in Fig. |l|, 

^measured _ q Hz and /,pi„ = 0.016 Hz). 

To quantify the presence of stripes in the maps of resid- 
uals, we can compute the value of the "striping" estimator 
rms(a^'^'^)/rms(a^^'^), because stripes tend to appear 
as structure grossly parallel to the iso- longitude circles. In 
the case of pure white noise with a uniform sky coverage, 
this value is 1. Here, because of the scanning strategy, the 
sky coverage is not uniform and the value of this estimator 
is greater than 1 , showing that it is not specific of the strip- 
ing. In order to get rid of the effect of non-uniform sky cov- 
erage, we express the estimator rms(a^'^'^)/rms(a^^'^) 
in units of rms(a^^' ' )/rms(a£Q ' ) for the white noise. 
This new estimator is specific to the striping. The results 
in Table ^ show the improvement achieved by the destrip- 
ing algorithm although the result is still not perfect. 



Fig. 11. The MoUweide projection of the residuals of the 
Q — Stokes parameter after destriping of 1// noise plus 
white noise with 77 = 1 

Fig. 12. The MoUweide projection of the residuals of the 
Q — Stokes parameter after zeroing the average of the cir- 
cles, for 1// noise plus white noise with 77 = 1. Although 
the remaining structures seem small, they are responsible 
for the excess of power in Cf , see Fig. [T5| . 

Table 1. Values of rms(a^'^''^)/rms(a^"^'"^) 
as a function of 77 = /knoo//spin in units of 
rms(a[^'^'^(WN))/rms(aJo^^-^(WN)) We have checked 
that the systematic difference between the zero-averaged 
E and B fields is randomly in favor of E and B depending 
on the particular sky simulation. 



Method 


/knee/ /spin 


T 


E 


B 


wliite noise 





1 


1 


1 


destriped 


0.5 


1.19 


1.05 


1.03 


zero-averaged 


0.5 


51.9 


9.23 


3.64 


undestriped 


0.5 


6.98 


7.04 


15.2 


destriped 


1 


1.24 


1.12 


1.19 


zero-averaged 


1 


52.2 


9.91 
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6. Discussions and Conclusions 

Comparison with other methods. Although no other 
method has yet been developped specifically for destriping 
polarized data, many methods exist for destriping unpo- 
larized CMB data, which could be adapted to polarized 
data as well. 

We first comment on the classical method which con- 
sists in modelling the measurement as 

mt = AtpTp + nt (24) 

where A is the so-called " pointing matrix" , T a vector 
of temperatures in pixels of the sky, mt the data and rit 
the noise. The problem is solved by inversion, yielding an 
estimator of the signal: 

fp = [A*N-^A]-^A^N'''^m, (25) 

where N = (nn*) is the noise correlation matrix and ^* is 
the transposed matrix of A. 
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Fig. 13. EfRciency of destriping for the T field with 
/knoc//spin ~ 1- The Sample variance associated to a pure 
white noise mission is plotted as the dotted lines. The 
"destriped spectrum" is very close to the white noise spec- 
trum (within the limits due to the sample variance). The 
zero-averaged and the "not destriped" spectra are a cou- 
ple of orders of magnitude above. The solid line represents 
a standard CDM temperature spectrum and the dashed 
line represents a CDM temperature spectrum with reion- 
ization. 
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Fig. 15. EfRciency of destriping for the E field for 
/knco//spin = 1. The zcro-avcragcd spectrum is not as bad 
as for T but the residual striping we can see in Fig. [l^ 
leads to some excess of power for low values of £ (up to 
£ ^ 100). We do not see such effect in the destriped spec- 
trum (and maps). The spectra for the B fields are very 
similar. The solid line represents a standard CDM E spec- 
trum and the dashed line represents a CDM E spectrum 
with reionization. 
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Fig. 14. EfRciency of destriping for the T field with 
/knco//spin = 5. The modelling of low- frequency noise with 
an offset is no longer sufficient and the destriping leaves 
some power at low values of £. Nevertheless, it remains 



a very i good way to significantly reduce the effect of low- 
freque ncy noise. 



This method can be extended straightforwardly to po- 
larized measurements, at the price of extending by a fac- 
tor of 3 X ft, the size of the matrix Atp, by 3 that of vec- 
tor Tp (replaced by {Ip,Qp,Up)), and by h that of the 
data stream (remember that h is the number of polarime- 
ters). The implementation of this formally simple solu- 
tion may turn into a formidable problem when megapixel 
maps are to be produced. Numerical methods have been 



Fig. 16. Same as Fig. ^but for /knco/ fspu 
tra for the B fields are very similar. 



5. The spec- 



proposed by a variety of authors (Wright, 1996; Tegmark 



1997), that use properties of the noise correlation matrix 



(symmetry, band-diagonality) and of the pointing matrix 
(sparseness) . Such methods, however, rely critically on the 
assumption that the noise is a Gaussian, stationary ran- 
dom process, which has been a reasonable assumption for 
CMB mission as COBE where the largest part of the un- 
certainty comes from detector noise, but is probably not so 
for sensitive missions such as Planck. Our method requires 
only inverting a 3n x 3n matrix where n is the number of 
circles involved, and does not assume anything on the sta- 
tistical properties of low frequency drifts. It just assumes 
a limit frequency (the knee frequency) above which the 
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Fig. 17. Same as Fig. |l^ for the £^T-correlation for 

f'knec/ /spin — 1- 



100.0 r 




Fig. 18. Same as Fig. |T| but for /kncc//spin = 5. 



smaller /knoc//spin the smaller the remaining correlation. 
Figs. 13-18 already contain the errors induced from the 
fact that we did not include these correlations in the co- 
variance matrix, and thus demonstrate that the effect is 

small for /knoo//spin ^ 1- 

Conclusion. The destriping as implemented in this paper 
removes low frequency drifts up to the white noise level 
provided that /knee //spin < 1- For larger /knee, the simple 
offset model for the averaged noise could be replaced with 
a more accurate higher-order model that destripes to bet- 
ter precision provided the scan s trategy allows to do so, as 
discussed in Delabrouille et al ( 1998a ). We are currently 
working on improving our algorithm to account for these 
effects. However, despite the shortcomings of our model, it 
still appears to be robust for small /knee and can serve as 
a first order analysis tool for real missions. In particular, 
our technique cannot only be used for the Planck HFI and 
LFI, but can also be adopted for other CMB missions with 
circular s canning strategies, such as COSMOSOMAS for 
instance ( Rebolo et al., 1999 ). 
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Appendix A: Destriping of unpolarized data 

We give here the formulas for the simpler case of destrip- 
ing temperature measurements with bolometers. The as- 
sumptions are the same as in the polarized case and we 
adopt the same notation for the common quantities. In- 
stead of polarimeters, we have h bolometers. Since the 
measurement is no longer dependent on the orientation of 
the bolometer, the model of the measurement is given by: 



M 



4 j'.t) 



noise can be considered as a white Gaussian random pro- 
cess. 

Another interesting method is the one that has been 
used by Ganga in the analysis of FIRS data ( |Ganga, 1994 ) , 
which is itself adapted from a method developed originally 
by Cottingham ( 1987 ). In that method, coefficients for 
splines fitting the low temperature drifts are obtained by 
minimising the dispersion of measurements on the pixels 
of the map. Such a method, very similar in spirit to ours, 
could be adapted to polarization. Splines are natural can- 
didates to replace our offsets in refined implementations 
of our algorithm. 

Here, we have assumed that the averaged noise can be 
modeled as circle offsets plus white noise (Eq. i.e. that 
the noise between different measurements from the same 
bolometer is uncorrelated after removal of the offset. This 
allowed us to simplify the to that shown in Eq. (^. 
In reality, the circular offsets do not completely remove 
the low frequency noise and there does remain some cor- 
relation between the measurements. The amount of cor- 
relation is directly related to the value of /kncc//spin; the 



Oi + /i.o s u with u 



(A.l) 



where Mi^^s is the /i— vector made of measurements by 
the h bolometers, Oi is an /i— vector containing the offsets 
for the i'th circle and hj^s is the temperature in the direc- 
tion of the intersection point labeled by {i,j,S}. u is an 
/i- vector, corresponding to the ^ matrix of the polarized 
case. The can be written as: 



i,iel(i),(5=±l 

N,-^ {M.^j^s - O, - h,,s u) . (A.2) 

In this case, the physical quantity uniquely defined at 
an intersection point is the temperature of the sky at 
this point. The constraint used here for removing low- 
frequency noise is then: 



(A.3) 
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Given this relation, the minimization of the with re- 
spect to Oi and li j s leads to the linear system: 



E 



X 1 I X ' 



(A. -A,) 



■ E 

iGl(i),5=±l 



where 

Xi — Ny^ u 



{\,,s - 'iJ,^,-s) (A.4) 



(A.5) 



corresponds to the Xi matrix of the polarized case, 

= -M^Ar,-iM„-5 (A.6) 

corresponds to the §ij-,5 local Stokes parameters of the 
polarized case and the scalar 



(A.7) 



corresponds to the 3-vector A^ of the polarized case. 

Temperature offsets appear through their differences 
Ai — Aj so this linear system is not invertible and we can 
use the same methods to invert the system than in the 
polarized case. The size of the matrix to invert is n x n. 
After evaluation of the offsets A,, one can recover the 
value of temperature for any sample k along circle i: 



li.k — 3j- fc — Ai 



Appendix B: Details of the simulations 



(A., 



The results presented in this paper correspond to the "Op- 
timized Configuration" involving 3 polarimeters and to 
an angular step of 18': the angle between two consecutive 
samples along a circle is 18' and there is one circle every 
18'. 

In this case, the number of circles is ric — 1200 and the 
number of samples on each circle is rig = 1195. The spin 
axis has a sinusoidal motion around the ecliptic plane with 
an amplitude of 8° and with 8 cycles during the mission. 
The aperture angle of the circles is 85°. The simulated 
noise has /knoc = f?/spin with -q = 1 and 5. 

The white noise variance is calculated based on the 
expected sensitivity on Q and U (3.7 fiK/K) of the (arbi- 
trarily selected) 143 GHz polarized channel of the Planck 
mission. 

All the maps (CMB, dipole, galaxy, simulation) are 
HEALPIX maps with 196 608 pixels of 27.48'. Only 12 
pixels (0.006% of the map) are not seen by the mission: 
their values are set to the average of the map. The sig- 
nal maps have been smoothed by a gaussian beam with a 
FWHM set to 2.5 x 18'. 

We have run other simulations with "Optimized Con- 
figuration" involving 3 or 4 polarimeters, and with cy- 
cloidal, sinusoidal or anti-solar spin axis trajectories (see 



Bersanelli et al ( |1996D and the Planck web page^ for 
additionnal information about proposed scanning strate- 
gies). The aperture angle of the circles have been taken 
in [70°, 75°, 80°, 85°, 90°]. In all these cases, the results 
are similar for /knee //spin ~ 1- For more pessimistic noise 
cases, the choice of the scanning strategy may have a 
strong impact on the quality of the final maps. A quan- 
titative study of this point is deffered to a forthcoming 
publication. 
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